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I.  INTRODUCTION 


Intermittently  since  World  War  I,  efforts  have  been  directed  toward 
the  research  and  development  of  a  device  capable  of  the  electromagnetic 
acceleration  of  macroparticles.  Such  a  device  has  the  potential  for  a 
number  of  military  applications,  probably  the  most  important  of  which  is 
the  so-called  electric  gun.  The  major  advantage  of  this  sort  of  gun 
over  the  more  conventional  ones  is  that,  at  least  in  principle,  much 
higher  muzzle  velocities  can  be  obtained.  In  addition,  the  guns  should 
be  relatively  smokeless  and  noiseless  and  their  projectiles  should 
experience  fairly  constant  acceleration. 

Although  the  efforts  mentioned  above  have  met  with  varying  degrees 
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of  success,  it  has  been  only  since  the  recent  experiments  at  the 
Australian  National  University  (ANU)  that  military  interest  in  the 
problem  has  been  rekindled.  These  experiments  demonstrated  conclusively 
that  the  500  MJ  homopolar  generator,  previously  developed  at  ANU,  could 
provide  sufficient  power  to  accelerate  a  mass  of  about  3g  to  a  velocity 
of  about  6  km/s  in  a  distance  of  about  3m.  The  renewed  military  inter¬ 
est  is  made  evident  by  the  current  DARPA/ARRADCOM  supported  project  under¬ 
way  at  Westinghouse  which  has  as  its  goal  the  acceleration  of  a  300g  mass 
to  a  velocity  of  3  km/s  in  a  distance  of  about  4m.  If  such  an  objective 
can  be  met,  the  development  of  an  electric  gun  may  become  a  feasible 
military  goal. 

A  number  of  devices  exist  for  producing  the  electromagnetic  accel- 
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eration  of  projectiles,  the  most  notable  being  the  rail  gun  and  the 
4  5 

mass  driver  9  .  The  experiments  at  ANU  and  those  now  underway  at 
Westinghouse  were  performed  using  a  rail  gun  and  it  is  with  this  device 
that  our  analysis  deals.  A  schematic  of  the  rail  gun  is  shown  in 
Figure  1.  Sides  one  and  two  represent  the  rails  which  carry  current 
from  the  source,  side  three,  to  and  from  side  four.  The  projectile  is 
represented  by  the  shaded  part  of  the  figure  and  may,  if  desired,  be  a 


1.  S.C.  Rashleigh  and  R.A.  Marshall,  Electromagnetic  Acceleration  of 
Macroparticles  to  High  Velocities’1 2 3 4 5,  J.  Appl.  Phys.  £9,  2540  (1978). 

2.  R.A.  Marshall,  ’’The  Australian  National  University  Rail  Gun 
Project”,  Atomic  Energy,  16,  January  1975. 

3.  J.P.  Barber,  ’’The  Acceleration  of  Macroparticles  and  a  Hypervelocity 
Electromagnetic  Accelerator”,  Ph.D.  Thesis  (Australian  National 
University,  1972)  (Unpublished). 

4.  H.  Kolm,  ’’Basic  Coaxial  Mass  Driver  Reference  Design”,  Third 
Princeton/AIAA  Conference  on  Space  Manufacturing  Facilities,  Paper 
No.  77-534,  Princeton,  NJ,  May  1977. 

5.  K.  Fine,  ’’Basic  Coaxial  Mass  Driver  Construction  and  Testing”, 

Third  Princeton/AIAA  Conference  on  Space  Manufacturing  Facilities, 
Paper  No.  77-535,  Princeton,  NJ,  May  1977. 
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Figure  1.  Model  for  rail  gun. 
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conductor  which  completes  the  circuit.  There  are  some  advantages,  how¬ 
ever,  to  having  the  projectile  an  insulator  and  allowing  the  current  to 
be  conducted  by  a  plasma  arc,  shown  by  the  region  between  x  =  xq  and 

x=£.  The  primary  advantages  are  that  better  contact  with  the  rails  can 
be  maintained  by  an  arc  than  by  a  solid  armature  and  that  the  stress 
imparted  to  the  projectile  is  likely  to  be  more  uniform.  In  practice 
the  arc  can  be  created  by,  for  example,  an  exploding  wire.  The  earliest 
use  of  the  arc-driven  rail  gun  was  apparently  in  experiments  by  Brast 

and  Sawle^  who  succeeded  in  accelerating  projectiles  of  a  few  mg  to 
velocities  of  a  few  km/s.  At  any  rate  the  current  distribution  pro¬ 
duces  a  magnetic  field  in  the  space  bounded  by  the  four  sides  which 
interacts  with  the  current  through  the  armature,  accelerating  it  in 
the  x  direction. 

3 

Some  theoretical  analysis  of  the  solid-armature  rail  gun  has  been 

7 

carried  out  and,  more  recently,  McNab  has  undertaken  the  more  difficult 
task  of  obtaining  some  estimates  of  the  properties  of  the  arc  in  the 
arc-driven  gun.  Although  these  calculations  provide  considerable  in¬ 
sight  into  the  problem  as  well  as  good  qualitative  results,  they  are 
based  upon  a  number  of  assumptions  which  are  unlikely  to  be  strictly 
valid  for  the  problem  at  hand.  It  is  assumed  in  McNab' s  analysis,  for 
example,  that  the  pressure,  density,  and  temperature  within  the  arc  are 
uniform  although  it  is  unlikely  that  all  three  flow  variables  will  be 
constant  in  an  arc  that  does  not  move  at  constant  velocity.  In  addition, 
it  is  assumed  that  atoms  which  constitute  the  arc  are,  at  most,  singly 
ionized.  However,  the  results  of  the  calculation  indicated  that  the 

o 

temperature  of  the  arc  was  about  57, 000° K,  and  one  can  show  that  at 
this  high  temperature  most  of  the  atoms  should  be  ionized  to  a  higher 
degree . 

The  inadequacy  of  these  assumptions  in  no  way  detracts  from  McNab' s 
work  since,  as  pointed  out,  the  results  do  appear  to  be  qualitatively 
reliable  and  the  calculations  are  rather  straightforward  to  carry  out. 

In  the  present  analysis,  however,  we  wish  to  extend  the  work.  In 
particular,  we  will  solve  the  fluid-dynamical  equations  in  the  accelerat¬ 
ing  arc  and,  thereby,  account  for  the  position  dependence  of  the  flow 
variables  where  necessary.  In  addition,  it  will  be  assumed  that  the 
plasma  may  consist  of  doubly  ionized  as  well  as  singly  ionized  atoms. 

The  results  will  be  applied  to  the  Rashleigh-Marshall  (RM)  experiment^ 
and  compared  briefly  to  the  results  obtained  by  McNab.  The  question  of 
how  the  arc  should  be  scaled  for  a  rail  gun  of  the  size  proposed  in  the 
Westinghouse  experiment  will  then  be  examined,  and  the  properties  of  the 
arc  under  those  conditions  will  be  investigated. 

6.  D.E.  Brast  and  D.R.  Sawle,  "Study  of  a  Rail-Type  MHD  Hypervelocity 
Projectile  Accelerator",  Proc.  Seventh  Hypervelocity  Impact 
Symposium,  1964,  Vol.  1,  p.  187  (unpublished). 

7.  I.  McNab,  "Electromagnetic  Acceleration  by  a  High  Pressure  Plasma" 

J.  Appl.  Phys.  5JU  2549  (1980'). 

8.  Authors'  unpublished  calculations. 
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The  report  is  organized  as  follows.  In  Sec.  II,  the  model  and 
assumptions  are  discussed  briefly.  In  Sec.  Ill,  the  formal  analysis 
of  the  problem  is  presented,  including  the  derivation  of  a  set  of 
coupled  equations  which  must  be  solved  to  determine  the  properties  of 
the  arc  and  the  motion  of  the  projectile.  In  Sec.  IV,  a  limiting-case 
analytic  solution  to  the  equations  is  found  and  their  numerical  solu¬ 
tion  for  the  more  general  case  is  discussed.  The  results  are  then 
applied  to  both  the  RM  experiment  and  to  the  proposed  Westinghouse 
experiment.  Finally,  in  Sec.  V,  some  additional  discussion  of  the  re¬ 
sults,  the  assumptions,  and  possible  future  investigations  is  given. 


II.  MODEL  AND  ASSUMPTIONS 

The  model  whose  properties  we  wish  to  consider  consists  of  an  arc- 
driven  rail  gun  such  as  shown  in  Figure  1.  The  rails  are  assumed  to  be 
perfectly  conducting,  infinitely  thin  sheets  which  are  infinitely  ex¬ 
tended  in  the  z  direction,  as  is  the  arc.  The  gun  itself  is  therefore 
two  dimensional,  although  the  fluid-mechanical  properties  of  the  arc 
are  assumed  to  vary  only  in  the  direction  of  propagation.  Certainly 
the  major  variation  does  occur  in  that  direction  because  of  the  purely 
mechanical  effect  of  accelerating  the  arc.  However,  in  a  real  system 
some  variations  occur  in  directions  normal  to  the  propagation  direction, 
and  these  are  not  accounted  for  in  the  model.  The  finite  conductivity 
a  of  the  arc  is  accounted  for,  and  the  assumption  that  the  rails  are 
perfect  conductors  is  justified  on  the  basis  that,  in  reality,  the 
conductivity  of  the  rails  is  several  orders  of  magnitude  higher  than 
that  of  the  arc.  The  potential  difference  V  which  is  applied  along  side 
three  of  the  model  is  assumed  to  vary  with  time  so  as  to  produce  a  con¬ 
stant  current  per  unit  height  j  on  the  surface  of  the  rails.  In  prac¬ 
tice,  such  a  condition  can  be  approximated  by  using  as  a  source  an  in¬ 
ductive  store  of  sufficiently  high  inductance. 


We  will  be  interested  in  solving  the  equations  which  govern  the 
properties  of  the  arc  and  the  acceleration  of  the  projectile  in  the 
steady-state  approximation.  By  steady  state  we  mean  that  the  accelera¬ 
tion  of  the  projectile  and  all  parts  of  the  arc  is  the  same  at  any  given 
time,  and  that  the  fluid-mechanical  properties  of  the  arc  are  independ¬ 
ent  of  time  in  a  frame  which  accelerates  with  the  arc.  It  is,  of 
course,  also  assumed  that  the  mass  of  the  arc  is  constant  in  time  so 
that  any  vaporization  of  the  rails  is  neglected  during  the  time  of 
acceleration.  A  number  of  specific  criteria  which  must  be  satisfied 
to  justify  the  assumption  of  steady  state  are  discussed  in  the  ensuing 
pages.  In  particular,  this  assumption,  as  well  as  the  others,  is  dis¬ 
cussed  in  some  detail  in  Sec.  V. 
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III.  ANALYSIS 


A.  Electrodynamics . 


The  first  problem  which  must  be  solved  is  to  determine  the  electro¬ 
magnetic  fields  associated  with  the  arc  gun  ^n  Figure  1.  From  the 
symmetry  of  the  problem  and  from  the  results  of  similar  calculations, 
we  assume  solutions  to  Maxwell's  equations  of  the  form 

£  =  E(x)  a 

.  /  (3.1) 

B  =  B(x)  a 

z 

where  1:  and  ^  are  the  electric  and  magnetic  induction  fields,  respec¬ 
tively.  The  appropriate  Maxwell  equations  then  become* * 


9_E  =  -3B 
3x  3t 


and 


1! 

3x 


=  -yJ 


(3.2a) 


(3.2b) 


where  y  is  the  magnetic  permeability  and  J  the  current  density  to  be 
determined  within  the  arc.  In  writing  Eq.  (3.2b)  we  have  neglected  the 
displacement  current  which  is  negligible  for  nonrelativistic  velocities 
of  the  arc  and  projectile. 

Equation  (3.2b)  implies  that  ?  is  constant  except  within  the  arc, 
and  since  it  must  vanish  at  infinity,  it  vanishes  everywhere  to  the 
right  of  x  -  £,  to  the  left  of  x  =  0,  and  above  and  below  the  perfectly 
conducting  rails.  Within  the  arc  Eq.  (3.2b)  implies 


B 


J  dx  +  c 


(3.3) 


where  c  is  an  integration  constant  which  can  be  determined  from  the 
continuity  of  B  and  from  the  boundary  condition 


BOO 


(3.4) 


We  have,  therefore, 


9.  R.M.  Fano,  L.J.  Chu,  and  R.B.  Adler,  Electromagnetic  Fields,  Energy, 
and  Forces  (Wiley,  New  York,  1960),  Chap.  9. 

*  Unless  otherwise  noted,  MKS  units  are  used  throughout  this  report. 
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J  dx, 


(3.5) 


B  = 


for  x  <  x  <  l  . 
o 

In  the  space  between  x=0  and  x^x^,  B  is  again  constant  and  to  satis¬ 
fy  the  boundary  condition  at  x=x  we  must  have 


rl 

B  =  y  J  J  dx  .  (3.6) 

x 

o 

However,  the  integral  on  the  right  of  Eq.  (3.6)  must  just  be  equal  to  j, 
the  current  per  unit  height  on  the  rails,  if  current  is  to  be  conserved 
in  the  steady  state.  Consequently,  we  have 


B  =  V3,  (3.7) 

for  0  <  x  <  x 

o 

The  electric  field  E  can  be  determined  most  simply  by  making  the 
change  of  variable 


x;  =  x  -  x  (t) 
o  v  J 


(3.8) 


in  Eq.  (3.2a).  For  steady-state  conditions  we  find 


where 


3E  _  3B 
9  x7-  "  V 


v 


dx 

o 

dt 


(3.9) 

(3.10) 


is  the  common  velocity  of  the  arc  and  projectile.  Integrating  Eq.  f3  9) 
we  find  ’ 


E  =  vB  +  F (t)  (3.11) 

where  F(t)  is  an  arbitrary  function  to  be  determined. 

The  function  F(t)  can  be  evaluated  by  noting  that  the  current  den- 
sity  i-n  the  arc  must  obey  the  constitutive  relation 

J  =  a(E  -  vB)  .  (3.12) 
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Thus,  F(t)  is  given  simply  by 

F (t)  =  J/ct  .  (3.13) 

Furthermore,  for  the  assumed  steady  conditions,  the  quantity  J/a  is 
independent  of  time  so 


E  =  Eq  +  vB,  (3.14) 

where  the  constant 

Eq  =  J/ct  (3.15) 

represents  the  electric  field  required  to  drive  the  current  through  the 
arc.  If  we  make  use  of  the  current-conservation  relation  mentioned 
earlier,  namely, 

i 

J(x)dx  =  j,  (3.16) 

o 


Eq  can  be  written  in  terms  of  known  quantities  as 


E 

o 


j 


o  (x)dx 


o 


(3.17) 


In  the  remaining  discussion  we  will  need  the  values  of  ^  and  only 
to  the  right  of  x=xq  and  between  the  perfectly  conducting  rails.  Writing 

these  in  terms  of  the  dimensionless  distance  £  defined  by 


5  =  (x  -  xo)/£a 


(3.18) 


we  have 

*  A 

B  =  y£E  j  ct(£)  d£  a  0  <  £  <  1 

aoj  z 

5  (3.19) 

$  =  0.  £  >  1  . 


In  these  equations  is  the  length  of  the  arc  measured  along  the  x 
axis,  viz.. 


£ 

a 


t-x  . 
o 


(3.20) 
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Similarly, 


£  =[  E0  +  vB(5)] 


(3.21) 


for  5  >  0  where  Eq  is  given  by  Eq.  (3.17). 

The  solutions  obtained  in  this  section  can  be  seen  to  satisfy  all 
boundary  conditions  appropriate  to  the  problem  and  the  assumption  made 

in  Eq.  (3.1)  is  therefore  justified.  A  more  arduous  derivation10,  empl 
ing  the  vector  potential  and  avoiding  the  assumption,  produces  the  same 
result. 

B.  Fluid  Mechanics. 

Having  obtained  the  values  of  the  electric  and  magnetic  fields,  we 
now  wish  to  derive  a  set  of  equations  appropriate  for  determining  the 
fluid-mechanical  properties  of  the  arc.  The  equations  which  express 
conservation  of  mass,  momentum,  and  energy  for  the  arc  as  a  whole  can 
1 1 

be  written 


f 


(3.22) 


where  p  is  the  mass  density  in  the  arc,  v  the  mean  flow  velocity,  e 
the  specific  internal  energy,  q  the  heat  flux,  f  the  force  per  unit 
volume  acting  on  the  plasma,  and  P  the  appropriate  element  of  the 

AA 

pressure  tensor  for  one-dimensional  considerations.  Equations  (3.22) 
are  identical  to  those  generally  encountered  in  fluid  mechanics  except 
for  the  force  f  which  is  electromagnetic  in  origin  and  except  for  the 

term  J  /a  which  represents  the  Joule  heating  of  the  plasma.  In  compli¬ 
ance  with  the  assumptions  previously  discussed,  the  flow  variables  have 
been  assumed  to  vary  only  in  the  direction  of  propagation  of  the  arc. 

It  is  clear  that  the  three  equations  in  Eqs.  (3.22)  cannot  uniquely 
specify  the  variables  P,  p,  v,  and  e  and  one  needs,  in  addition,  an 
equation  of  state.  We  assume,  therefore,  that  the  plasma  obeys  an 
ideal-gas  law,  namely, 

TOT  Authors'  unpublished  calculations. 

11.  A.B.  Cambel,  Plasma  Physics  and  Magnetof luidmechanics  (McGraw-Hill 
New  York,  1963),  Chap.  8.  "  “ 
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(3.23) 


V 

P  =  (1  +cO  P  — 
o 

where  kD  is  Boltzmann's  constant,  T  the  temperature,  m  the  atomic  mass 
B  u 

o£  the  ion  or  neutral,  and  a  the  ratio  of  electrons  to  heavy  particles 
in  the  plasma.  The  method  for  determining  a  will  be  discussed  in  the 
following  section.  All  flow  variables  are  assumed  to  vary  only  slowly 
with  position  so  that  the  gas  is  in  a  state  of  local  thermal  equili¬ 
brium. 

Further  analysis  of  Eqs.  (3.22)  is  most  easily  accomplished  in  a 
frame  moving  with  the  velocity  v  of  the  plasma  arc.  We  therefore  make 
the  change  of  variable  defined  in  Eq.  (3.18).  Noting  that  in  the  moving 
frame  v  is  independent  of  £  and  that  p,  e,  and  P  are  independent  of  time, 
we  find  that  the  first  of  Eqs.  (3.22)  is  satisfied  identically,  while  the 
remaining  two  become 

p**r-!£=f  (3-24) 

a 

and 

|f  =  l  J 2/a  .  (3.25) 

a 

In  Eq.  (3.24)  cl  represents  the  common  acceleration  of  the  arc  and  pro¬ 
jectile  and  P  has  been  replaced  by  the  pressure,  P,  since  no  velocity 

gradients  are  assumed  to  exist. 

The  force  per  unit  volume  f  acting  on  the  plasma  arises  from  the 
interaction  of  the  magnetic  induction  field  within  the  arc  and  the 
current  through  it.  The  force  acts  in  the  positive  x  direction  and 

its  magnitude  is  given  by  ^ 


f  =  JB  . 


Consequently,  Eq.  (3.24)  can  be  written  in  the  form 

£ 


P(0  =  K 


J(E)  BOO  -  p(0«  I  ds  +  c 


(3.26) 


(3.27) 


where  c  is  an  integration  constant  which  can  be  determined  from  the 
condition  that  the  pressure  at  the  back  of  the  projectile  must  be  suf¬ 
ficient  to  provide  it  with  an  acceleration  CL.  Specifically,  if  we  let 


12.  J.D.  Jackson,  Classical  Electrodynamics  (Wiley,  New  York,  1962), 
Chap .  5 . 
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Pz  represent  the  density  of  the  projectile  multiplied  by  its  length 
P 

along  the  x  axis,  we  must  have 

1 

P(5=l)  =  J  (JB-pa)d£  +  c  =  p  a  .  (3.28) 

o  p 


It  is  clear,  however,  that  the  term  JB,  integrated  over  the  entire  vol¬ 
ume  of  the  arc,  yields  the  total  force  acting  on  the  arc  and  projectile 


(p 


a 

P 


)a  = 


i 

a 


JCO  B(5)  d£ 


(3.29) 


where  represents  the  average  arc  density  multiplied  by  its  length 

a 

along  the  x  axis.  In  the  appendix  it  is  proved  that  the  acceleration  a 
is  independent  of  how  the  current  is  distributed  in  the  arc;  in  fact, 
the  simple  expression 


a  = 


(3.30) 


is 

at 


obtained.  At 
last 


P«)  = 


any  rate,  Eqs.  (3.28)  and  (3.29)  imply  c=0,  so  we  have 


l 

a 


[J  (£)  B(0  -  P(5)a]  dg. 


(3.31) 


The  second  term  on  the  right-hand  side  of  Eq.  (3.31)  is  clearly  negligi¬ 
ble  with  respect  to  the  first  provided  <<  as  it  is  under  most 


conditions . 

From  Eq.  (3,23)  the  mass  density  of  the  arc  is  given  by 


m  P 
o 


K  kgT(l  +  a) 

and  the  length  of  the  arc  can  be  determined  from  the  condition 

,1 


(3.32) 


./ 


P (5)d£  =  p. 


(3. 33) 
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Solving  for  Z  , 

3. 


we  obtain 


(.3.34) 


The  parameter  p  is,  of  course,  independent  of  Z  . 

x,  a 

a 

Provided  the  temperature  of  the  arc  is  sufficiently  high,  heat  is 
transferred  within  the  arc  and  to  its  surroundings  primarily  by  radia¬ 
tion  and  conduction  can  be  neglected.  Therefore,  q  in  Eq.  (3.25)  can 
be  taken  to  represent  the  heat  flux  due  to  radiation,  and  we  have  upon 
integration 


(3.35) 


The  constant  of  integration  c  can  be  determined  from  the  fact  that  at 
the  interface  between  the  arc  and  the  (assumed)  vacuum  behind  it  (5=0), 

13 

the  flux  must  approximately  reduce  to 

q (0)  =  -  2asTo4,  (3.36) 


where  Tq  is  the  temperature  of  the  arc  at  5  =  0  and  cs  is  Stefan's  con- 

-  8  J 

stant,  namely,  5.67  X  10  - -t—  .  Consequently,  we  have 

m  deg  s 


q  = 


2  4 

J  ,  d£  -  2o  T  , 
a  so 


(3.37) 


where  the  unknown  temperature  Tq  must  be  determined. 

If  the  radiative  mean  free  path  within  the  arc  is  small  compared 
to  the  arc  length,  as  we  shall  find  to  be  the  case,  the  flux  q  can  be 

13 

approximated  by  the  "radiative  heat  conduction"  equation. 


13  Y.B.  Zel'dovich  and  Y.P.  Raizer,  Physics  of  Shock  Waves  and  High- 
Temperature  Hydrodynamic  Phenomena  (Academic,  New  York,  1966), 

Vol .  I,  Chap.  2. 
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(3.38) 


3  l  d£ 

3. 


where  X  represents  the  mean  free  path  for  radiation  averaged  over  the 
entire  frequency  spectrum.  Obviously  to  solve  the  equation  an  additional 
boundary  condition  must  be  specified.  To  obtain  the  condition  we  note 
that  during  the  entire  time  of  acceleration  energy  is  transferred  to  the 
projectile  surface,  continuously  heating  its  interior.  As  the  projectile 
temperature  rises,  the  transfer  of  energy  from  arc  to  projectile  becomes 
less  efficient,  and  this  change  in  the  transfer  rate  with  time  will 
clearly  produce  a  time-dependent  effect  on  the  temperature  in  the  arc. 
Thus,  in  a  strict  sense,  the  temperature  profile  of  the  arc  is  not  com¬ 
pletely  steady.  It  is  clear,  however,  that  the  projectile  temperature 
at  any  time  is  substantially  less  than  that  of  the  arc  and,  consequently, 
little  error  should  be  made  in  assuming  that  the  arc  at  the  projectile 
surface  is  maintained  at  some  constant  temperature.  Furthermore,  since 
the  radiative  mean  free  path  is  small  compared  to  the  arc  length,  the 
properties  of  the  arc,  except  very  near  the  boundary,  should  be  inde¬ 
pendent  of  this  temperature.  For  these  reasons  it  is  convenient  to 
assume  that  the  projectile  absorbs  all  the  radiation  which  is  incident 
upon  it,  just  as  does  the  vacuum  behind  the  arc,  and  therefore  the  flux 
at  £=1  reduces  to 


q(l)  =  2  asT14 


(3.39) 


where  Tj  is  the  arc  temperature  at  the  projectile  surface.  This  con¬ 
dition  supplies,  then,  the  additional  boundary  condition.  Further 
discussion  of  this  assumption  is  given  in  Sec.  V. 

We  have,  upon  integrating  Eq.  (3.38), 


(3.40) 


From  Eqs .  (3.37),  (3.39),  and  (3.40),  we  obtain  for  the  unknown  tempera¬ 
tures  Tq  and  T  the  relationships 


(3.41) 


o 


and 


n 


(3.42) 
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Solving  for  Tq  and  T  and  substituting  in  Eqs.  (3.37)  and  (3.40) 
we  find 


q  = 


l 

a 


and 


T 


(3.43) 


(3.44) 


Equations  (3.43)  and  (3.44)  provide  two  coupled  equations  for  the 
heat  flux  q  and  temperature  T.  Formally  the  equations  can  be  solved  by 
the  transformation 


q  =  A  +  F(S)  . 


(3.45) 


After  substitution  into  Eq.  (3.43)  we  find 


FC5)  =  *a 


(3.  46) 


where  <  1/A  >  is  the  inverse  of  the  mean  free  path  averaged  over  the 
length  of  the  arc,  i.e., 


< 


(3.48) 


Equation  (3.44)  can  now  be  written 
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C.  Degree  of  Ionization. 

The  parameter  a,  occurring  in  Eq.  (3.23),  must  be  obtained  before 
a  solution  to  the  fluid-mechanical  equations  considered  in  the  last 
section  can  be  found.  In  practice,  this  parameter  results  from  the 
solution  of  a  hierarchy  of  equations  generally  referred  to  simply  as 
the  Saha  equation.  In  this  section  we  outline  the  method  by  which  the 
solution  is  obtained.  For  greater  detail  the  reader  is  referred  to  the 
14 

literature 

We  recall  that  the  parameter  a  represents  the  ratio  of  the  number  of 
electrons  to  the  number  of  heavy  particles,  N„.  We  let  x.  be  the  ratio 

H  3 

of  the  number  of  atoms  ionized  j  times  (referred  to  as  "jth  ions")  to 
Nr*  IT  is  then  evident  that  a  is  given  by 


a  =  l  j  x.  (3.50) 

j  J 

and  the  electron  number  density  in  the  arc  by 


n 

e 


pot 

m 

o 


(3.51) 


The  Xj's  are  obviously  normalized  such  that 

I  X.  =  1.  (3.52) 

3 

In  Eqs.  (3.50)  and  (3.52),  the  sums  over  j  run  over  all  values  for  which 
appreciable  ionization  occurs. 


It  has  been  shown  that  the  x^'s  satisfy  the  following  system  of  equa¬ 
tions,  usually  called  the  Saha  equation: 


x.  .a 
3  +  l 

x. (1+a) 
3 


3/2 

(kgT) 5/2 


/  kgT 


K.+1(T,P).(3.53) 


In  Eqs.  (3.53)  h  is  Planck's  constant  divided  by  2tt , 
Ij  The  ionization  energy  needed  to  ionize  the  atom  j 
are  electronic  partition  functions  for  the  jth  ion. 
given  by 


mg  the  electron  mass, 

times,  and  the  Z . 's 
3 

Specifically,  Z.  is 


14.  See  Ref.  13,  Chapter  3. 
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Vc)  -l  *ji e 


•Un/kBT(S) 


(3.54) 


where  IL  ^  is  the  energy  of  the  electronic  state  of  the  ion  and 

where  is  the  appropriate  degeneracy  factor  for  this  level.  The 

partition  functions  depend  upon  position  in  the  arc  since  the  tempera¬ 
ture  is  a  function  of  position  in  the  problem  under  consideration. 
Similarly,  the  x^'s  are  also  position  dependent. 

Equations  (3.53)  are  deceptively  simple  looking.  In  practice  the 
equation  for  any  always  depends  on  xj+^  and,  thus,  the  system  of 

equations  must  be  terminated  at  some  point.  If,  for  example,  we  assume 
that  only  single  ionization  occurs,  then. 


a  =  x^,  (3.55) 

and  only  the  term  corresponding  to  j=0  is  retained  in  Eqs.  (3.53).  The 
single  resulting  equation  becomes 


1  -x 


~2  =  K^T.P) 


(3.56) 


if  we  use  Eqs.  (3.50)  and  (3.52).  If  desired,  the  equation  can  be  solved 
for  x^  (or  a),  and  substituted  into  the  fluid-dynamic  equations  of  the 

previous  section. 


For  the  problem  under  study  here,  however,  it  is  usually  insuffic¬ 
ient  to  assume  only  single  ionization.  In  the  event  that  double  ioni¬ 
zation  is  also  allowed. 


a  =  X1  +  2x2 

and  Eqs.  (3.53)  become 

x:(x  +2x2) 

(1-x1-x2)(1+x1+2x2)  =  K1(T’P) 

x  (x  +2x  ) 

Vl'+x^x.)  °  K2(T-P>- 


(3.57) 


(3.58) 
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Again  we  have  used  Eqs.  (3.50)  and  (3.52).  These  equations  are  most 
easily  solved  by  iteration.  After  considerable  algebra,  they  can  be 
rewritten  in  the  form 


X1  = 


K1X2 

~!kZ 


1 

4K2(1-x) 

1  -L  - 

1/2  J 
-  1 

K1X2 

\ 

L 

) 

(3.59) 


K1 

4K(l-x)' 
1  +  - 

1/2 

1  K,X- 

1  2 

(2+3K1)x2 

1+K, 


4K1(l+K1)(l+x2-2xp 


x2(2+3K1)2 


+1 


1/2 


(3.60) 


A  value  for  x^  is  then  assumed,  substituted  into  the  right-hand  side  of 
Eq.  (3.60)  and  a  new  value  of  x^  calculated.  The  process  is  then  re¬ 
peated  until  successive  values  agree  to  within  any  desired  accuracy.  The 
value  for  x1  then  follows  directly  from  Eq.  (3.59). 

To  illustrate  the  results  we  have  plotted  in  Figure  2  values  of 

xl  and  x2  as  a  function  of  log  T  for  a  pressure  of  10  N/m  .  The  plasma 

was  assumed  to  be  a  gas  of  partially  ionized  copper  atoms.  The  partition 

functions  in  Eq.  (3.54)  were  calculated  using  values^  for  the  g.^  and 

Uj  shown  in  Table  I.  Calculations  were  undertaken  retaining  first  ten, 

and  then  fifteen,  terms  in  the  partition  functions  and  no  discernible 
differences  were  found  in  the  results.  The  initial  value  of  x2  was 

chosen  to  be  0.01  and  successive  values  were  found  to  agree  to  within 

one  part  in  10  after,  at  most,  a  few  tens  of  iterations.  As  can  be 
seen  from  the  graph  no  appreciable  ionization  exists  below  log  T  *  9 
(T=8, 100°K) .  Second  ionization  is  negligible  below  log  T  «  10  (T  = 
22,000°K)  and  then  rises  fairly  rapidly.  Once  x2  has  become  large  (above, 

say,  log  T  *=  10.5  or  T  =  36,000°K),  triple  ionization  is  probably  also 
important,  but  is  not  accounted  for  in  the  analysis  above. 


D.  Summary  of  Governing  Equations. 

In  the  analysis  of  the  preceding  sections  we  have  derived  a  system 
of  coupled  equations  sufficient  to  determine  the  electromagnetic  fields, 
the  acceleration  of  the  projectile,  and  the  fluid-mechanical  properties 

15.  Data  taken  from  C.E.  Moore,  "Atomic  Energy  Levels",  National  Bureau 
of  Standards  Circ.  No.  467,  V.  II,  Washington,  DC,  1952. 
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log  T 


Figure  2.  Ionization  factors  as  a  function  of  temperature. 
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Atomic  Energy  Levels  and  Degeneracy  Factors  for  Copper  and 
Its  Ions.  Energies  are  expressed  in  °K  (from  Ref.  15). 
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of  the  arc.  For  easy  reference  the  most  pertinent  equations  are  listed 
below  [see  Eqs.  (3.15),  (3.17),  (3.19),  (3.30), •  (3.31),  (3.32),  (3.34), 
(3.45),  (3.46),  (3.47),  (3.49),  (3.50),  (3.51),  (3.59),  and  (3.60)]: 


J  =  a  E 

o 


E 

o 


m  P 
o 

P  "  kBT(l+x1+2x2) 


(3.61) 

(3.62) 

(3.63) 

(3.64) 

(3.65) 

(3.66) 


n 


e 


p (x1+2x2) 


m 

o 


(3.67) 

(3.68) 

(3.69) 

(3.70) 
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(3.71) 


K1X2 

2K_ 


1  + 


4K2(1-x2) 


K1X2 


1/2 


X2 


1  + 


4K2(1-x2) 

K1X2 


T7T 


-  l 


(2+3K1)x2 


1+K, 


4K1(l+K1)(l+x2-2xp 


X2  ^2+^K1) 2 


+  1 


1/2 


(3.72) 


-1 


In  writing  Eqs.  (3.61)  -  (3.72)  we  have  assumed  that  the  plasma  is 
at  most  doubly  ionized.  If  we  are  provided  with  appropriate  expressions 
for  the  conductivity  a  and  the  mean  free  path  X,  the  equations  are  suf¬ 
ficient  to  determine  the  twelve  unknowns:  J,  B,  P,  p,  T,  q,  n  ,  x  ,  x?, 

EqJ  a,  and  £,  .  The  first  nine  of  the  unknowns  are  position  dependent, 

whereas  the  last  three  are  constant. 


IV.  SOLUTION  OF  EQUATIONS 
A.  Limiting-Case  Analytic  Solution. 


In  this  section  we  will  be  concerned  primarily  with  the  numerical 
solution  of  the  equations  summarized  in  Sec.  HID  and  their  application 
to  the  experimental  data.  Before  turning  to  the  numerical  results, 
however,  we  will  present  a  limiting-case  analytic  solution  which  is 
valid  whenever  (1)  the  mass  of  the  arc  is  negligibly  small  with  respect 
to  the  mass  of  the  projectile,  (2)  the  conductivity  of  the  arc  is  con¬ 
stant,  and  (3)  the  radiative  mean  free  path  is  inversely  proportional 
to  the  pressure.  The  analytic  solution  serves  as  a  useful  check  for 
the  computer  program  used  to  evaluate  the  numerical  results  and  provides 
some  physical  insight  into  the  behavior  of  the  equations  as  well.  In 
addition,  the  above  assumptions  are  not  totally  unreasonable  for  the 
problems  which  will  be  considered  as  will  be  seen  when  the  more  general 
numerical  results  are  presented. 

For  constant  conductivity,  Eqs.  (3.61)  and  (3.62)  imply 

J-f~  C4-1) 

a 

so  the  current  is  uniformly  distributed  throughout  the  arc  as  expected. 
Equation  (3.63)  then  becomes 
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B  =  yj  (1-5) 


(4.2) 


and  the  magnetic  induction  field  varies  linearly  from  the  value  yj  at 
the  trailing  edge  of  the  arc  to  zero  at  the  leading  edge.  Similarly, 
if  we  neglect  the  second  term  on  the  right-hand  side  in  Eq.  (3.65)  in 
compliance  with  the  first  assumption  above  we  have 

P  =  yj2(5-?2/2).  (4.3) 

Under  assumption  (3)  above,  we  have  from  Eq.  (4.3) 

X 

X  =  - f—  (4.4) 

K-rl  2 

where  Xq  is  a  constant  independent  of  position  within  the  arc*.  Equa¬ 
tions  (3.46)  and  (3.47)  can  now  be  evaluated  and  yield  for  the  limiting- 
case  solution 


and 


F(C) 


2 

1 


i  a 
a  o 


(4.5) 


A  = 


8a  ft 
o  a 


(4.6) 


where  oq  denotes  the  constant  conductivity  of  the  arc.  In  obtaining 

Eq.  (4.6)  we  have  noted  that  <  1/X  >  is  simply  1/3Xq  when  Eq.  (4.4)  is 

applicable.  Provided  the  mean  free  path  is  small  compared  to  the  arc 
length,  Eq.  (4.6)  can  be  expanded  in  powers  of  X  /£  to  obtain 

O  cl 


A 


8a  ft 
o  a 


(4.7) 


*It  may  be  seen  that  X  »  as  £  ->  0,  contradicting  our  statement  that 
X  <<  £  within  the  arc.  As  we  shall  see,  however,  reasonable  values  of 

-4 

X  are  of  the  order  of  10  m  while  £  is  of  the  order  0.1m,  so  the  con- 
o  a 

dition  is  violated  only  very  near  the  boundary.  In  fact,  we  found  X<£ 

a 

over  98%  of  the  arc  in  our  numerical  calculations.  The  small  error 
incurred  in  the  calculation  resulting  from  the  lack  of  validity,  say, 
of  radiative  heat  conduction  in  this  small  region  is  deemed  negligible. 
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The  expressions  for  F(£)  and  A  can  now  be  used  to  determine  the 
heat  flux  q  and  temperature  T  within  the  arc.  The  flux  follows  directly 
from  Eq.  (3.68)  and  yields 

.2 

« =  iV  (?  -  5/8)  +  °W  (4-8) 

a  o 


where  °(X0Aa)  denotes  higher-order  terms  in  the  parameter  X  /£  .  Thus, 

O  3. 

the  flux  changes  sign  approximately  five-eighths  of  the  distance  between 
the  trailing  edge  of  the  arc,  where  it  is  negative,  and  the  leading  edge, 
where  it  is  positive.  The  asymmetric  nature  of  the  flux  results  from  the 
fact  that  the  mean  free  path  is  a  function  of  position. 


Substituting  Eqs.  (4.4)  and  (4.5)  into  Eq.  (3.69)  and  evaluating  the 
results  we  find 


T 


4 


.2 

J 

64a  a  X 

S  0  o 


20  jS.  +  15(1-4X  /5l aK2 
a 


-  21(1-4X ft/2Ha)53  +  6  ?4 

o  a 


°(x  V*  J 


(4.9) 


In  the  interior  part  of  the  arc  where  0  <  £  <  1,  the  terms  of  order 
X0Aa  can  be  neglected  in  Eq.  (4.9)  and  we  have 


T  = 


3j2g2(5-2Q(l-Q 

64X  a  a 
o  s  o 


(4.10) 


Here  the  temperature  is  independent  of  as  should  be  expected  for 

Xo^a  <<:  1 *  Very  near  the  Pr°j ectile,  however,  as  C  +  1,  the  lowest- 
order  terms  in  Eq.  (4.9)  vanish  and  terms  proportional  to  X  /£  must 
be  retained.  We  obtain  then 


161  a  o 
a  s  o 

Similarly,  as  E,  ->  0,  or  at  the  trailing  edge  of  the  arc,  we  find 
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(4.12) 


T  = 


II 


.2 


16£  a  o 
a  s  o 


Very  near  the  boundaries  of  the  arc,  then,  the  temperature  does  depend 
on  i  and  is  lower  than  that  at  the  center  of  the  arc  by  a  factor  of 

a  t- 

the  order  of  (X  /£  ) 4.  The  physical  meaning  of  these  results  will  be 

O  3. 

discussed  in  some  greater  detail  when  the  numerical  solutions  are  pre¬ 
sented. 


The  solutions  presented  above  should  offer  some  insight  into  the 
physical  meaning  of  the  results  when  we  attempt  to  interpret  the  numeri¬ 
cal  calculations.  Further  simplification  of  the  equations  in  Sec.  HID 
is  not  possible  without  additional  assumptions  and,  thus,  numerical 
techniques  must  still  be  employed  to  obtain  a  complete  solution  even  in 
the  simplified  case.  However,  if  one  can  make  a  further  assumption  that 
the  ion  concentrations  and  x9  are  constant  across  the  arc,  a  complete 

uncoupling  of  the  equations  is  possible.  Such  an  approximation  might 
be  reasonable,  for  instance,  in  a  temperature  range  in  which  either  single 
or  double  ionization  was  essentially  complete,  the  other  being  negligible. 
Under  this  assumption  we  have  from  Eqs.  (3.66)  and  (3.67) 

1 

i  m 

=  ° _ 

a  p£  kB(1+Xl+2x2} 
a 

o 

Since  Eq.  (4.10)  gives  the  temperature  in  the  arc  everywhere  except  for  a 
negligible  contribution  near  the  boundaries,  this  expression  may  be 
used  for  T  in  Eq.  (4.13)  and  we  have,  after  employing  Eq.  (4.3)  as 
well , 


/ 


P(Q 

T(5) 


dK 


(4.13) 


.1 

£ 


.3/2 
V)  m 


p£  kB(1+Xl+2x2} 


d?  . 


(4.14) 


The  integral  on  the  right-hand  side  of  Eq.  (4.14)  is  a  constant  inde¬ 
pendent  of  any  properties  of  the  arc  and  its  value  is  about  0.46. 

The  value  of  i  for  any  particular  set  of  input  parameters  can  be 
a 

calculated  from  Eq.  (4.14)  and  results  substituted  into  all  equations 

in  this  section  which  depend  on  £  .  Consequently,  under  these  admittedly 

a 

restrictive  conditions,  explicit  values  for  all  unknowns  discussed  in 
Sec.  IIID  can  be  obtained. 
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B.  Numerical  Solution. 


The  equations  in  Sec.  HID  have  been  written  in  a  form  which  makes 
them  particularly  amenable  to  solution  by  iterative  techniques.  Before 
proceeding,  however,  expressions  for  the  conductivity  a  and  radiation 
mean  free  path  X  are  required.  For  the  conductivity  we  have  used  a 

1  z:  _  -I  n 

result  obtained  by  Spitzer  and  coworkers  for  a  completely  ionized 

gas,  namely. 


<KO 


2.63  X  10"2y  T3//2 
E 

z 


Here  Z,  given  by 


log 


23X107T5/2 

Z  /~n~ 


Z 


V4x2 

V2x2 


(4.15) 


(4.16) 


is  a  position-dependent  number,  lying  between  1  and  2,  which  gives  an 
indication  of  the  degree  of  ionization  in  the  gas.  The  multiplicative 
factor,  Yg,  actually  depends  weakly  on  Z.  For  the  problem  at  hand, 

however,  the  gas  is  nearly  completely  doubly  ionized  and  y  can  be  well 

E 

approximated  by  the  value  0.6833. 


An  expression  for  the  radiative  mean  free  path  averaged  over  the 
frequency  spectrum  (Rosseland  mean  free  path)  has  been  derived  by 

19 

Raizer  for  a  multiply  ionized  plasma.  The  derivation  is  based  on 
the  assumption  that  the  atoms  and  ions  which  constitute  the  gas  can  be 
approximated  by  hydrogenic  atoms  or  ions  with  an  appropriate  charge; 
it  accounts  for  both  bound-free  and  free-free  (bremsstrahlung)  transi¬ 
tions.  One  finds 


X(5)  = 


0.91  X  10UkBT3(l+x1  +  2x2) 


-yv  -w  -w 

(l-x1-x2)e  +4xxe  “  B  +9x0e  3  B 


(4.17) 


for  the  case  in  which  the  gas  is  at  most  doubly  ionized. 


16.  L.  Spitzer,  Physics  of  Fully  Ionized  Gases  (Interscience,  New  York 

1965),  Chap.~5^  " 

17.  R.S.  Cohen,  L.  Spitzer,  and  P.  McR.Routly,  "The  Electrical  Conduc¬ 
tivity  of  an  Ionized  Gas",  Phys.  Rev.  80,  230  (1950). 

18.  L.  Spitzer  and  R.  Harm,  "Transport  Phenomena  in  a  Completely  Ionized 
Gas",  Phys.  Rev.  S9_,  977  (1953). 

19.  Y.P.  Raizer,  "Simple  Method  for  Computing  the  Mean  Range  of  Radia¬ 
tion  in  Ionized  Gases  at  High  Temperatures,"  Sov.  Phys.  -  JETP 

37,  769  (1960).  See  also  Ref.  13,  Chap.  5. 
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The  numerical  solution  of  the  equations  in  Sec.  HID  was  then 
carried  out  as  follows.  The  interval  0  <  5  <  1  was  divided  into  a  few 
hundred  equally  spaced  segments,  and  values  for  the  position-dependent 
unknowns  were  estimated  at  every  grid  point  from  the  approximate  solu¬ 
tion  in  Sec.  IVA.  Values  for  the  constants  £  and  E  were  also  esti- 

a  o 

mated  in  this  manner  and  a  was  calculated  from  Eq.  (3.30).  The  esti¬ 
mated  values  were  then  used  as  initial  values  in  an  iterative  technique 
by  substituting  them  on  the  right-hand  side  of  Eqs.  (3.61)  -  (3.72) 
and  calculating  new  values  on  the  left-hand  side.  All  integrals  were 
evaluated  simply  by  the  trapezoidal  rule.  For  the  second  iteration,  a 
weighted  average  of  the  old  and  new  values  was  used  on  the  right-hand 
side.  The  process  was  then  repeated  and  successive  values  were  found 

4 

to  differ  by  only  about  one  part  in  10  after  a  few  hundred  iterations. 
Calculations  were  performed  using  different  step  sizes  and  it  was  found 
that  approximately  400  grid  points  gave  sufficiently  accurate  results. 

To  the  extent  possible  the  mechanics  of  the  program  was  checked  using 
the  analytic  solution  in  Sec.  IVA. 

C .  Application  to  Rashleigh-Marshall  Experiment. 

As  pointed  out  in  Sec.  I,  Rashleigh  and  Marshall'*'  have  used  an  arc- 
driven  rail  gun  to  accelerate  a  3g  mass  to  a  velocity  of  about  6  km/s. 

We  now  wish  to  solve  the  equations  in  Sec.  IIID  using  data  appropriate 
to  the  RM  experiment  and,  thereby,  to  analyze  the  arc.  With  two  excep¬ 
tions,  the  data  which  are  listed  in  Table  II  are  available  either  di¬ 
rectly  from  the  experiment  or  can  be  easily  calculated  therefrom.  It 
has  been  necessary,  however,  to  estimate  the  arc  mass  since  this  quan¬ 
tity  was  not  measured  experimentally  and  we  have  chosen  to  use  the 

7 

value  of  O.lg  estimated  by  McNab  in  his  earlier  calculations.  We  have 
in  subsequent  work  varied  m  ,  holding  the  other  parameters  fixed,  and 

these  results  will  be  discussed  in  Sec.  IVD.  In  addition,  in  the  ex¬ 
periment  the  outer -most  edges  of  the  rails  were  surrounded  by  insulat¬ 
ing  material  which  excluded  the  arc  from  this  region.  Consequently, 
the  rail  height,  h  ,  was  somewhat  larger  than  the  arc  height  measured 

along  the  rails,  h^.  In  the  model  studied  here,  both  these  heights 

were  assumed  equal  and  we  have  therefore  used  for  h  an  effective  height 
found  from  the  geometric  mean  of  the  two  values.  The  remaining  quan¬ 
tities  in  the  table  are  largely  self-explanatory.  The  potential  Vq  is 

the  voltage  measured  across  the  muzzle  of  the  gun;  it  gives  rise  to 
the  electric  field  in  Eq.  (3,21)  and,  since  B=0  at  the  muzzle,  it 
represents  physically  a  purely  resistive  drop  across  the  plasma.  The 
atomic  mass  mQ  was  easily  found  to  be  the  value  indicated  for  the  cop¬ 
per  arc  used  in  the  experiment. 
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TABLE  II.  Experimental  Data  for  Rashleigh-Marshall  Experiment.  The 
data,  taken  mostly  from  Ref.  1,  were  used  in  the  numerical 
solution  of  the  governing  equations  in  Sec.  IIID. 


Quantity 

Description 

Value 

w 

Rail  separation 

1.27  X  lo"2m 

h 

P 

Plasma  height  on  rails 

E 

CM 

1 

O 

H 

X 

CM 

H 

h 

r 

Rail  height 

1.91  X  lo"2m 

h 

Effective  rail  height 

1.56  X  I0"2m 

m 

P 

Projectile  mass 

3  X  10" 3  kg 

i 

Pulsed  current 

3  X  105  A 

j 

Current  per  unit  height 

7 

on  rails 

1.92  X  10  A/m 

V 

0 

Muzzle  voltage 

160  V 

* 

m 

a 

Arc  mass 

10"4  kg 

m 

0 

Ion  (or  atom)  mass 

1.1  X  io~25  kg 

a 

Average  acceleration 

6  X  10^  m/s^ 

p£ 

a 

[See  Eq.  (3.29)] 

0.51  kg/m2 

P£ 

[See  Eq.  (3.29)] 

15.1  kg/m2 

P 

*  Estimated  value. 


The  equations  in  Sec.  IIID  have  been  solved  in  the  manner  discussed 
in  Sec.  IVB  and  using  the  data  in  Table  II.  Results  for  the  position- 
independent  unknowns  are  indicated  in  Table  III,  while  graphs  of  posi¬ 
tion-dependent  quantities  are  shown  by  the  solid- line  curves  in  Figures 
3-10.  (The  curves  denoted  by  xfs  will  be  discussed  in  Sec.  IVE.)  These 
variables  are  plotted  as  a  function  of  the  dimensionless  distance 
The  brackets  on  quantities  in  Table  III  indicate  that  these  variables 
were  averaged  over  the  length  of  the  arc. 

The  value  obtained  for  the  acceleration  is  about  a  factor  of  two 
higher  than  the  experimental  value  in  Table  II.  Since  a  is  a  sensitive 
function  of  j,  however,  close  agreement  cannot  be  expected  in  view  of 
the  rather  arbitrary  manner  in  which  h  was  chosen.  The  length  of  the 
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TABLE  III.  Results  of  Numerical  Calculation  for  RM  Experiment.  The 
arc  mass,  m  ,  was  O.lg. 

3. 


Quantity 


a 


Z 


a 


E 

o 

<  T  > 


<  Z> 


V 

o 

<  n  > 
e 


Value 

1.48  X  107  m/s2 
9.20  cm 

3. 72  X  103  volts/m 
5.61  X  104  °K 
1.93 

47  volts 
9.85  X  1025  m"3 


arc  Z^  was  not  explicitly  measured  in  the  experiment,  but  the  calcu¬ 
lated  value  is  consistent  with  the  estimate  of  around  10  cm. 

The  mean  temperature  is  probably  somewhat  high,  and  it  is  likely 
that  the  high  value  results  from  the  one-dimensional  character  of  the 
model  in  which  energy  has  been  allowed  to  radiate  only  from  the  ends  of 
the  arc,  not  from  the  sides.  In  a  two-or-three-dimensional  model  more 
surface  area  is  available  for  radiation  and  the  temperature  must  fall. 

A  crude  estimate  of  the  mean  temperature  in  a  three-dimensional  model 
may  be  obtained  as  follows.  In  one  dimension,  the  surface  area  avail¬ 
able  for  radiation  is 


A^=2hw  , 


(4.18) 


while  in  three  dimensions  one  has 

A  =  2  hw  +  2  H  w  +  2  £  h  (4.19) 

j  a  a 

Using  data  for  the  RM  experiment  and  assuming  negligible  change  in  Z  , 

a 

we  find  A^/A^  «  14.  Since  for  a  blackbody,  the  temperature  varies 

inversely  as  the  fourth  root  of  the  available  surface  area,  we  might 
expect  a  drop  in  <  T >  of  about  a  factor  of  two. 
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The  average  of  the  parameter  Z,  defined  in  Eq.  (4.16),  is  included 
in  the  table  in  order  to  give  some  indication  of  the  degree  of  ioniza¬ 
tion  in  the  arc.  Evidently,  most  of  the  ions  are  doubly  ionized  with 
very  few  singly  ionized  atoms  in  most  of  the  arc.  The  high  degree  of 
second  ionization  makes  it  likely  that  some  triple  ionization,  not 
accounted  for,  probably  also  exists.  Again,  however,  Z  should  be  con¬ 
siderably  less  in  a  two-  or  three-dimensional  model  because  of  the 
lower  temperature. 

Finally,  the  muzzle  voltage  which  is  easily  calculated  from  the  re¬ 
lation 

Vo  =  Eo  w  (4.20) 

was  found  to  be  47  volts,  about  a  factor  of  three  smaller  than  that  mea¬ 
sured  in  the  experiment.  At  the  lower  temperatures  expected  in  three 
dimensions,  Vq  should  increase  since,  according  to  Eq.  (4.15),  the  con¬ 
ductivity  of  the  arc  decreases  with  decreasing  T. 

We  see  therefore  that  the  calculated  quantities  which  can  be  com¬ 
pared  to  the  experimental  results  agree  to  within  about  a  factor  of 
three.  The  agreement,  we  feel,  is  reasonable  in  view  of  the  restrictive 
assumptions  made  in  the  model  as  well  as  the  arbitrary  way  in  which  m 

and  h  were  obtained. 

Plotted  in  Figure  3  is  the  pressure  as  a  function  of  position  in 
the  arc.  The  pressure  rises  from  zero  at  the  trailing  edge  of  the  arc 
(a  possible  definition  of  the  trailing  edge)  to  a  value  at  the  projec¬ 
tile  sufficient  to  provide  an  acceleration  a,  namely,  a.  For  the 

P 

appropriate  data,  this  value  is  roughly  220  MPa,  or  about  2200  atmos¬ 
pheres.  If  the  current  density  in  the  arc  were  uniform,  P  would  vary 
2 

as  / 2  (see  Sec.  IVA) .  The  nonconstant  current  density,  however, 

gives  rise  to  deviations  from  this  type  of  behavior  observed  in  the 
graph.  Similarly,  the  magnetic  induction  field,  B,  shown  in  Figure  4, 

varies  from  about  24  T  at  the  trailing  edge  to  zero  at  the  projectile. 

The  variation  would  be  linear  for  constant  current  density  as  can  be  seen 
from  Eq.  (4.2). 

The  asymmetric  nature  of  the  temperature  profile,  shown  in  Figure 
5,  results  from  the  position  dependence  of  the  radiation  mean  free  path. 
Specifically ,  toward  the  leading  edge  of  the  arc  where  the  density  in¬ 
creases,  the  mean  free  path  decreases.  Consequently,  photons  created 
m  the  right -most  part  of  the  arc  are  radiated  away  with  more  difficulty 
than  those  in  the  left-most  part,  and  there  is  a  high  concentration  of 
energy  near  the  projectile.  Very  near  the  projectile  the  temperature 
drops  suddenly  as  was  predicted  by  the  analytic  results  in  Sec.  IVA. 
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DISTANCE 

Figure  3.  Pressure,  in  MPa,  as  a  function  of  position  in  the  arc: 

- ,  calculation  in  Sec.  IVC;  xxx,  calculation  in 

Sec.  IVE. 
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DISTANCE 


Figure  4.  Magnetic  induction  field,  in  Tesla,  as  a  function  of 

position  in  the  arc:  - ,  calculation  in  Sec.  IVC;  xxx, 

calculation  in  Sec.  IVE. 
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Figui e  5.  Temperature,  in  thousands  of  degrees  Kelvin,  as  a  function  of 

position  in  the  arc:  - ,  calculation  in  Sec.  IVC;  xxx, 

calculation  in  Sec.  IVE. 
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Clearly  once  the  photons  are  within  a  mean  free  path  or  so  of  the  pro¬ 
jectile  surface,  they  are  able  to  "see”  the  boundary  and  are  radiated 
outward  with  little  difficulty.  Consequently,  the  temperature  drops 
dramatically  over  this  very  short  distance.  It  is  also  interesting  to 

note  that  the  temperature  at  £  =  0,  2.53  X  io4  °K,  is  larger  than  that 

at  the  projectile,  2.16  X  io4  °K,  by  a  factor  of  about  1.17.  This  fac¬ 
tor  is  very  close  to  that  predicted  by  Eqs.  (4.11)  and  (4.12),  namely, 

(5/3) 1/4  *  1.14. 

The  mass  density,  p,  shown  in  Figure  6,  varies  nearly  linearly  from 

<7 

the  expected  value  of  zero  at  the  trailing  edge  to  about  10  kg/m  (about 
eight  times  the  density  of  air  at  STP)  at  the  leading  edge.  Very  near 
the  boundary  a  sudden  rise  in  the  density  by  about  a  factor  of  five  is 
observed.  The  rise  results  because  as  the  temperature  drops  at  the 
boundary,  the  degree  of  ionization  also  drops,  and  both  these  effects 
tend  to  reduce  the  pressure  near  the  projectile.  Thus,  a  steady  pres¬ 
sure  at  the  projectile  can  be  maintained  only  if  there  is  a  proportionate 
increase  in  mass  density  at  the  surface. 

In  Figure  7  is  plotted  the  electron  density,  ng,  which  varies  in 

much  the  same  manner  as  the  mass  density  as  should  be  expected.  The 
rise  near  the  projectile  surface,  however,  is  not  so  rapid  as  for  the 
mass  density;  even  though  there  are  more  ions  near  the  boundary,  they 
are  less  highly  ionized  than  those  interior  to  the  arc  because  of  the 
lower  temperature. 

The  curve  for  the  current  density,  J,  shown  in  Figure  8,  is  largely 
self-explanatory.  The  current  falls  to  zero  at  the  trailing  edge  of 
the  arc  where  the  electron  density  is  zero,  increases  with  increasing 
temperature  across  the  arc  as  should  be  expected  from  Eq.  (4.15),  and 
finally  drops  as  T  drops  near  the  boundary.  It  should  perhaps  be  noted 

from  Eq.  (4.15)  that  a  varies  nearly  as  T^^,  depending  only  weakly  on 
the  remaining  parameters.  (Z  is  nearly  constant  over  most  of  the  arc 
as  can  be  seen  in  Figure  9.) 

In  Figure  9  is  shown  the  mean  ionic  charge,  Z,  defined  in  Eq.  (4.16). 
Except  very  close  to  the  projectile,  we  found  that  essentially  all  the 
atoms  in  the  arc  -were  at  least  singly  ionized  with  very  few  neutrals 
present.  Therefore,  within  the  arc  a  value  of  Z  =  2  corresponds  to  x^^  = 

0  and  x^  =  1 ,  while  a  value  of  Z  =  1.8  corresponds  roughly  to  x^  =  0.33, 

X2  =  0.67.  The  unusual  behavior  of  Z  near  the  trailing  edge  of  the  arc 

results  from  the  factor  1/P  in  the  Saha  equation.  This  dependence  re¬ 
flects  the  fact  that  electron  reattachment  is  a  two-body  interaction, 
whereas  ionization  does  not  necessarily  require  a  collision.  However, 
the  results  are  dependent  on  the  assumption  of  thermal  equilibrium  which 
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Figure  6.  Mass  density,  in  kg/m  ,  as  a  function  of  position  in  the  arc: 

- ,  calculation  in  Sec.  IVC;  xxx,  calculation  in  Sec. 

IVE. 
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l igure  7.  Electron  density,  in  m  and  normalized  by  the  constant  fac- 
25 

tor  10  ,  as  a  function  of  position  in  the  arc:  - , 

calculation  in  Sec.  IVC;  xxx,  calculation  in  Sec.  IVE . 
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Figure  8.  Current  density,  in  MA/m  ,  as  a  function  of  position  in  the 

arc:  - ,  calculation  in  Sec.  IVC;  xxx,  calculation  in 

Sec.  IVH. 
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Figure  9.  Mean  ionic  charge  as  a  function  of  position  in  the  arc: 

- ,  calculation  in  Sec.  IVC;  xxx,  calculation  in  Sec. 

IVE. 
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is  not  valid  near  5=0,  owing  to  the  small  density  there.  Therefore, 
the  detailed  behavior  of  Z  (and  the  other  parameters)  very  near  5=0  is 
probably  open  to  question.  In  any  case,  the  region  in  doubt  is  only  a 
small  portion  of  the  overall  profile  and  should  have  little  effect  on 
the  mean  properties  and  overall  dynamics  of  the  arc  as  suggested  pre¬ 
viously.  Slightly  to  the  right  of  5=0,  the  degree  of  ionization  is 
determined  predominantly  by  the  temperature  and  increases  with  increasing 
temperature  across  the  arc.  Very  near  the  projectile,  the  temperature 
drops  dramatically  with  only  a  small  variation  in  the  pressure,  so  Z 
likewise  drops.  In  fact,  at  the  projectile  surface,  x^  =  0.38  and  x0  ** 

0  so  only  38%  of  the  heavy  particles  which  constitute  the  arc  are  at 
all  ionized. 

Finally,  we  show  in  Figure  10  the  radiation  heat  flux  q.  For  the 
limiting-case  analytic  solution,  we  found  [see  Eq.  (4.8)],  a  linear  be¬ 
havior  in  q  as  a  function  of  5  with  the  flux  changing  sign  at  5  **  5/8 

**  0.63.  For  the  more  general  calculation,  the  behavior  is  fairly  well 
approximated  by  a  linear  variation  and  the  sign  changes  around  5=0.72. 

D.  Dependence  of  Results  on  Arc  Mass. 

Since  it  was  necessary  to  estimate  the  value  of  the  arc  mass  in 
the  previous  section,  it  is  of  some  interest  to  ask  how  the  results  de¬ 
pend  cn  m  .  We  have  therefore  performed  the  calculation  for  a  number 
a 

of  different  values  of  m  ,  holding  the  remaining  parameters  fixed  at  the 
values  in  the  previous  section. 

Typical  of  the  results  are  those  shown  in  Table  IV  which  are  for 
an  arc  mass,  m  ,  of  0.05  g.  As  can  be  seen  from  the  table,  the  accelera- 

cl 

tion  of  the  arc  is  very  slightly  higher  than  before  because  of  the  de- 


TABLE  IV.  Variation  of  RM  Results  With  Arc  Mass.  Results  were  obtained 
using  imput  data  in  Table  II,  except  that  m  was  given  by 
0 . 05g .  a 


Quantity 

Value 

a 

1.51 

x  io7 

m/s2 

l 

a 

4.54 

cm 

E 

0 

7.47 

X  io3 

volts/m 

<  T  > 

5.66 

X  io4 

°K 

Z 

1.94 

V 

o 

95  volts 

<  n  > 

1.00 

y  ,26  -3 

X  10  m 

e 
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Figure  10. 


GJ 

Heat  flux,  in  - -  ,  as  a  function  of  position  in  the  arc: 

s  m 

- ,  calculation  in  Sec.  IVC;  xxx,  calculation  in  Sec. 

IVE. 
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crease  in  the  mass  of  the  arc.  Primarily  the  acceleration  is  determined 
by  the  more  massive  projectile,  however,  and  the  arc  mass  has  only  a 
minimal  effect.  The  length  of  the  arc  has  decreased  by  approximately  a 
factor  of  two  as  might  be  expected  intuitively.  With  the  decrease  in 
l  ,  Eq  and  consequently  Vq  rise  by  about  a  factor  of  two  resulting 

from  the  increase  in  resistance  of  the  arc.  Likewise,  the  current  den¬ 
sity  J  in  the  arc,  though  not  shown  in  the  table,  is  everywhere  higher 
by  about  a  factor  of  two,  owing  to  the  decreased  area  of  the  arc. 

It  is  noteworthy  that,  within  the  interior  of  the  arc,  the  remain¬ 
ing  pcsition-dependent  variables  considered  previously  (B,P,T,p  .n^x^ 

and  x,.)  depend  only  very  weakly  on  m  when  plotted  as  a  function  of  the 

dimensionless  distance  £.  In  fact,  when  m  was  changed  from  0.1  g  to 

a 

0.05  g  we  observed  changes  of,  at  most,  a  few  percent  across  the  inter¬ 
ior  of  the  arc,  though  more  significant  changes  did  occur  very  near  the 
boundaries.  Except  for  the  temperature,  this  result  is  perhaps  not  too 
surprising.  It  does  appear,  however,  that  the  higher  rate  of  energy 
dissipation  in  the  smaller  arc  would  surely  result  in  a  higher  tempera¬ 
ture.  The  reason  for  the  unexpected  result  is  that,  even  though  the 
Joule  heating  is  indeed  greater  in  the  smaller  arc,  so  are  the  tempera¬ 
ture  gradients  across  the  arc.  Therefore,  although  energy  is  dissipated 
at  a  higher  rate,  it  is  also  radiated  away  more  easily.  For  the  one¬ 
dimensional  model  under  study  here,  the  two  effects  almost  exactly  com¬ 
pensate  each  other  so  the  arc  seeks  nearly  the  same  steady-state  tempera¬ 
ture  in  the  two  cases.  This,  incidentally,  will  probably  not  be  the  case 
in  a  two-  or  three-dimensional  model  and  one  expects  the  temperature  to 
rise  with  decreasing  m  (or  Z  ) . 

cl  SL 

The  manner  in  which  the  parameters  discussed  above  vary  with  can 

be  seen  approximately  from  the  special-case  solution  in  Sec.  IVA.  Thus, 
one  sees  from  Eqs.  (4.2)  and  (4.3)  that  B  and  P  do  not  vary  with  m&  as 

a  fur.ction  of  £ .  From  Eq.  (4.14),  on  the  other  hand,  Z  varies  directly 

a 

with  m  (or  p  )  in  this  limiting  case  and,  therefore,  as  m  decreases  by 
a  *a 

a  factor  of  two,  so  does  Z  .  In  accordance  with  Eq.  (4.1),  then,  J 

a 

should  vary  inversely  with  m  as  should  E  and  V  .  T,  on  the  other 
J  a  oo 

hand,  is  nearly  independent  of  Z  or  m  within  the  interior  of  the  arc 

a  a 

[Eq.  (4.10)]  but  some  variation  can  be  expected  very  close  to  the 
boundaries  [Eqs.  (4.11)  and  (4.12)].  Finally,  since  p,  Xy  and  x2  de¬ 
pend  only  on  P  and  T,  one  expects  them  also  to  vary  significantly  with 
m&  only  very  near  the  boundaries. 

Of  course,  the  variations  predicted  by  the  limiting-case  solution 
are  only  approximate  in  view  of  the  assumptions  made  in  obtaining  that 
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solution.  It  is  interesting,  nevertheless,  that  the  approximate  scaling 
laws  predicted  from  results  in  Sec.  IVA  appear  to  hold  very  nearly  even 
for  solutions  to  the  more  general  equations. 

E •  Application  to  the  Proposed  Westinghouse  Experiment  and 

Scaling  Factors. 

We  have  also  undertaken  calculating  the  properties  of  the  arc  in  a 
gun  of  the  size  proposed  for  the  Westinghouse  experiment.  In  particular, 
the  input  data  for  that  calculation  are  shown  in  Table  V.  The  rail  gap, 
rail  height,  and  projectile  mass  were  chosen  to  be  values  appropriate 
for  the  proposed  experiment;  the  value  of  the  pulsed  current  i  was  chosen 
so  as  to  provide  an  acceleration  sufficient  to  reach  a  velocity  of  3 
km/s  in  a  distance  of  4m;  and  the  arc  mass  was  selected,  mostly  by  trial 
and  error,  to  give  an  arc  length  of  the  order  of  10  cm. 


TABLE  V.  Input  Data  for  Large  Gun.  These  data  were  used  in  the 

numerical  solution  of  the  governing  equations  in  Sec.  HID; 
they  correspond  roughly  to  a  gun  comparable  in  size  to  that 
proposed  in  the  Westinghouse  experiment 


Quantity 

Description 

Value 

w 

Rail  separation 

5  X  10'2  m 

h 

Rail  height 

5  X  10-2  m 

m 

P 

Projectile  mass 

0.3  kg 

i 

Pulsed  current 

7.35  X  105  A 

j 

Current  per  unit  height  on  rails 

1.47  X  107  A/m 

m 

a 

Arc  mass 

8  X  10-4  kg 

p* 

a 

[See  Eq.  [3.29)] 

0.32  kg/m2 

p* 

p 

[See  Eq.  [3.29)] 

120  kg/m2 

General  results  of  the  numerical  calculation  are  shown  in  Table  VI 
and,  as  can  be  seen,  quantities  are  of  the  same  order  of  magnitude  as 
obtained  in  the  RM  calculation.  The  largest  difference  in  the  two  sets  of 
data  occurs  for  the  potential  across  the  rails  which  is  larger  than  that 
for  the  RM  experiment  (with  ma  =  O.lg)  by  nearly  a  factor  of  five.  The 
higher  value  in  the  present  case  results  from  the  higher  current  and 
higher  arc  resistance  in  the  larger  gun. 


48 


Table  VI.  Results  of  Numerical  Solution  for  Large  Gun.  Input  data 
used  in  the  calculation  are  given  in  Table  V. 


Quantity  Value 


a 

1.13  X  106 

m/s2 

l 

a 

7.82  cm 

E 

0 

4. 14  X  103 

volts/m 

<  T  > 

4.68  X  104 

°K 

<  Z  > 

1.89 

V 

0 

207  volts 

<  n  > 

25  -  3 

7.11  X  10  m  J 

e 

Graphs  of  position-dependent  quantities  are  denoted  for  this  calcula¬ 
tion  by  the  curves  drawn  with  x!s  in  Figures  3-10.  The  general  discus¬ 
sion  of  the  curves  for  the  RM  experiment  applies  here  as  well,  and  no 
further  discussion  will  be  given. 

It  is  of  interest  to  ask  whether  we  can  determine  some  approximate 
scaling  laws  from  the  analytic  solution  in  Sec.  IVA  and,  if  so,  to  what 
extent  these  laws  predict  the  results  obtained  for  the  larger  gun. 

Fairly  simple  scaling  laws  can  be  derived  from  the  analytic  solution 
provided  we  make  three  additional  assumptions:  First,  we  assume  that 
the  degree  of  ionization  does  not  appreciably  change  from  one  experi¬ 
ment  to  the  other  as  can  be  seen  to  be  the  case  for  the  RM  and  Westing- 
house  experiments.  Second,  we  neglect  the  variation  of  the  logarithmic 

3/2 

term  in  Eq.  (4.15),  i.e.,  assume  a  T  .  Finally,  we  note  that,  for 

constant  ionization,  X  in  Eq.  (4.17)  varies  approximately  as  T^'  ^/p^.  ^ 
If  we  make  use  of  these  assumptions,  some  simple  algebra  yields  the 
scaling  laws  shown  in  Table  VII.  The  pertinent  equations  for  deriving 
the  appropriate  scaling  factor  are  indicated  in  the  right-most  column 
of  the  table.  To  use  the  table,  the  value  of  a  particular  quantity  in 
the  RM  experiment  (column  1)  is  multiplied  by  the  ratio  of  scaling  fac¬ 
tors  for  the  large  gun  and  RM  gun.  The  numerical  value  of  this  ratio 
for  the  two  experiments  is  shown  in  column  3.  The  result  is  then  the 
predicted  value  of  the  quantity  for  the  larger  gun. 
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Table  VII.  Approximate  Scaling  Factors.  Scaling  factors  were  derived 
from  the  limiting-case  analytic  solution  in  Sec.  IVA. 


Quantity 

Scaling  Factor 

Value 

Equations 

£ 

a 

, .16/11 

pe  h 

0.92 

3.66,4.3,4.10, 

a 

4.14,  4.15 

a 

j2/(P*  +  p£  ) 

0.076 

3.64 

a  p 

V 

o 

.18/11  , 

1  w/p 

4.08 

3.62,3.66,4.3, 

a 

4.10,4.14,4.15, 

4.20 

E 

0 

.18/11 

J  /pS> 

1.03 

3.62,3.66,4.3, 

a 

4. 10,4.14,4.15 

B 

j 

0.77 

4.2 

P 

.2 

J 

0.59 

4.3 

J 

.27/11  . 

3  /p£ 

0.84 

3.66,  4. 1,4. 3, 

a 

4.10,4.14,4.15 

T 

j  6/11 

0.87 

3.66,4.3,4.10, 

4.15 

P 

j 16/11 

0.68 

3.66,4.3,4.10, 

4. 15 

n 

e 

j 16/11 

0.68 

3.66,  3.70,  4.3, 
4.10,4.15 

q 

)29/n/„t 

0.80 

3.66,4.3,4.8, 

a 

4.10,4.14,4.15 
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Values  obtained  in  this  manner  for  the  parameters  in  Table  VI  were: 

a  =  1.13  X  106  m/s2,  £a  =  8.46  cm,  Eq  =  3.83  X  103  volts/m,  <  T  >  = 

4.88  X  104  °K,  V  =  192  volts,  and  <  n  >  =  6.70  X  1025m"3.  As  can  be 

o  e 

seen,  values  of  the  parameters  obtained  from  the  scaling  laws  agree  with 
the  actual  theoretical  values  in  Table  VI  to  within,  at  worst,  about  15%. 
We  have  also  investigated  the  validity  of  the  scaling  laws  in  predicting 
the  values  of  position-dependent  quantities  and,  again,  found  reasonable 
agreement.  Typical  of  the  results  is  that  shown  in  Figure  11.  The  curve 
denoted  by  the  xfs  is  the  theoretical  value  of  the  arc  temperature  for 

the  larger  gun  plotted  as  a  function  of  £  (see  Figure  5);  the  curve  de¬ 

noted  by  the  0!s  is  the  same  function  but  obtained  from  the  scaling  laws. 
The  latter  curve  was  deduced  by  multiplying  the  temperature  at  every 
point  £  in  the  RM  experiment  by  the  appropriate  scaling  factor,  namely, 
0.87.  As  can  be  seen  the  results  agree  everywhere  to  within  less  than 
10%. 


It  appears,  therefore,  that  no  major  difficulties  should  be  en¬ 
countered  if  it  should  be  desirable  to  employ  an  arc  in  the  proposed 
Westinghouse  experiment.  All  properties  of  the  arc  are  expected  to  be 
the  same  order  of  magnitude  as  for  the  RM  experiment.  Furthermore, 
approximate  values  of  the  pertinent  parameters  can  be  obtained  from  the 
scaling  factors  in  Table  VII. 


V.  DISCUSSION 

In  this  final  section,  we  will  discuss  the  calculation  briefly, 
particularly  the  major  assumptions,  indicate  the  limitations  of  the 
model,  and  suggest  what  future  calculations  might  be  desirable.  To 
summarize,  we  have  proposed  a  model  for  determining  the  steady  accelera¬ 
tion  of  the  arc  and  projectile  in  the  arc-driven  electric  gun  and  for 
describing  the  fluid-mechanical  properties  of  the  arc.  A  set  of  twelve 
coupled  equations  has  been  derived  which,  when  solved,  yields  these 
properties.  The  theory  has  been  applied  to  the  RM  experiment  as  well 
as  to  the  experiment  proposed  by  Westinghouse.  It  has  been  demonstrated 
that  use  of  an  arc  in  the  larger  gun  should  pose  no  major  obstacles. 

As  pointed  out  in  Sec.  I,  some  earlier  analysis  of  arc  dynamics  in 

7 

the  RM  experiment  has  been  carried  out  by  McNab  .  He  assumed  that  one- 
third  of  the  measured  muzzle  voltage  53  volts)  occurred  across  the 
plasma  and  used  this  experimental  value  to  calculate  the  temperature  of 
the  arc.  The  spatial  variation  of  all  flow  parameters  was  neglected 
and  the  gas  was  assumed  to  be,  at  most,  singly  ionized.  In  addition, 
the  experimental  value  of  the  acceleration  was  used  in  the  calculation 
and  the  arc  length  was  assumed  to  be  10  cm.  Values  obtained  for  the 
pressure,  temperature,  electron  density,  and  arc  mass  were:  P=110  MPa, 
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Figure  11.  Comparison  of  numerical  results  of  temperature  profile  with 
results  obtained  from  scaling  laws:  xxx,  numerical  results 
from  calculation  in  Sec.  IVE;  000,  results  obtained  from 
scaling  laws  in  Table  VII. 
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T  =  5.7  X  104deg,  ne  =  7  X  1025m  3,  and  ma  =  O.lg.  Clearly,  McNab's 

results  are  in  order-of-magnitude  agreement  with  results  obtained  in 
Sec.  I VC. 

The  rationale  for  assuming  that  only  one-third  of  the  potential 
difference  between  the  rails  occurs  across  the  arc  is,  as  McNab  pointed 
out,  that  some  contact  potential  surely  exists  at  the  rail-arc  inter¬ 
face.  This  contact  potential  has  not  been  included  in  the  present  cal¬ 
culation  primarily  because  it  is  not  clear  experimentally  how  signifi¬ 
cant  a  potential  drop  actually  exists.  Admittedly,  however,  the  con¬ 
tact  potential  is  probably  not  negligible  and  some  further  investigation 
is  appropriate.  Presumably,  the  effect  could  be  included  in  the  model 
with  little  additional  difficulty. 

The  major  assumption  that  has  been  employed  in  the  present  calcu¬ 
lations  is  that  the  entire  acceleration  process  can  be  approximated 
by  conditions  appropriate  to  the  steady  state.  This  approximation 
should  be  reasonable  provided  the  time  required  for  a  steady  state  to 
be  established  is  small  compared  to  the  total  time  of  acceleration. 
Obviously,  the  validity  of  the  approximation  cannot  be  rigorously  ascer¬ 
tained  without,  in  fact,  solving  the  time-dependent  problem.  Neverthe¬ 
less,  several  necessary  conditions  might  be  mentioned. 

First,  it  is  clear  that  for  the  steady-state  approximation  to  be 
valid,  the  current,  which  initially  flows  on  the  surface  of  the  arc, 
must  diffuse  through  it  in  a  time  small  compared  to  the  total  accelera¬ 
tion  time.  Alternatively,  the  skin  depth,  6,  of  the  plasma  must  be 
large  compared  with  its  length.  The  skin  depth  of  the  plasma  at  time 

20 

t  can  be  estimated  from  the  expression 


where  oq  here  denotes  the  mean  conductivity  of  the  arc.  The  second 

equality  in  Eq.  (5.1)  follows  from  Eqs.  (3.15)  and  (4.20)  with  J  ap¬ 
proximated  by  j/l  .  For  an  acceleration  time  appropriate  to  the  RM 

experiment,  namely,  1.6  ms  and  assuming  Vq  =  47  volts,  the  theoretical 

value,  we  find,  6  **  30  cm.  Clearly,  5  is  considerably  larger  than  the 
arc  length  of  9.2  cm.  That  the  skin  depth  here  is  larger  than  those 
generally  associated  with  solid  conductors  results  from  the  lower  con¬ 
ductivity  in  the  gaseous  state. 


20.  W.R.  Smythe,  in  American  Institute  of  Physics  Handbook,  edited 
by  D.E.  Gray  (McGraw-Hill,  New  York,  1957),  p.  5-85. 
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A  second  condition  which  must  be  satisfied  in  order  to  justify  the 
steady-state  approximation  is  that  the  time  required  for  the  plasma  to 
reach  its  steady-state  internal  energy  must  be  small  compared  to  the 
total  acceleration  time.  The  internal  energy  of  the  plasma  at  any  time 
can  be  determined  from  the  relation 


Eint 


ET  +  E.  +  E 
Tie 


(5.2j 


where  E^,  E^  and  Eg  denote  the  thermal,  ionization,  and  electronic 

energies,  respectively,  whereas  the  energy  which  has  been  dissipated  in 
the  arc  at  time  t  is  given  by 

'  P2 

(5.3) 
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Assuming  that  the  gas  is  completely  doubly  ionized,  a  simple  estimate 
reveals  that  for  the  RM  experiment,  E.nt  *  6  kJ.  If  we  neglect  losses 

in  energy  due  to  radiation  during  the  time  for  which  the  plasma  is  being 
heated,  the  time  required  for  the  plasma  to  reach  its  steady  internal 
energy  can  be  found  by  substituting  E.^  on  the  left-hand  side  of  Eq. 

(5.3)  and  solving  for  t.  We  find  t  »  0.4  ms,  a  value  somewhat  smaller 
tain  the  total  acceleration  time  of  1.6  ms.  The  assumption  that  radia- 
tion  losses  are  negligible  during  the  time  of  heating  is  probably  rea¬ 


sonable,  since  the  radiation  flux  varies  as  T4.  It  is  of  interest  to 
note,  incidentally,  that  the  total  energy  dissipated  in  the  arc  during  the 
acceleration  time  of  1.6  ms  is  about  20  kJ.  This  energy  is  considerably 
less  than  the  kinetic  energy  of  the  arc  and  projectile  which  is,  using 
the  theoretical  value  of  the  acceleration,  about  870  kJ.  In  the  model, 
the  remaining  energy  supplied  by  the  source  is  stored  in  the  magnetic 
field,  and  it  can  be  shown  that  this  energy  is  identical  to  the  kinetic 
energy  of  the  arc  and  projectile. 


Finally,  a  less  stringent  condition  which  must  be  satisfied  in  order 
to  justify  the  steady-state  approximation  is  that  the  time  of  accelera- 
ion  must  be  large  compared  to  the  time  necessary  for  a  sound  wave  to 
traverse  the  arc.  If  this  condition  is  not  met,  it  is  clear  that  the 
flow  variables  cannot  become  equilibrated  during  the  time  of  accelera¬ 
tion,  and  that  they  will  be  time  dependent  in  a  frame  accelerating  with 
t  ie  arc.  One  can  estimate,  however,  that  the  sound  speed  in  the  plasma 
is  of  the  order  of  a  few  kilometers  per  second.  Thus,  the  time  necessary 
or  a  disturbance  to  traverse  an  arc  10  cm  long  is  probably  somewhat  less 

than  10  s.  Again,  this  is  substantially  less  than  the  acceleration  time 
of  1.6  ms . 


The  above  considerations  do  not  prove  the  validity  of  the  steady 
assumption  and  further  investigation  is  desirable.  In  principle,  a 
time-dependent  theoretical  treatment,  accounting  for  the  the  initial 
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compression  and  expansion  of  the  arc  and  the  diffusion  of  current  within 
it,  should  be  possible.  Such  a  calculation,  however,  would  be  immensely 
more  difficult  than  that  undertaken  here  and  would  necessitate  additional 
experimental  data  and  assumptions.  It  would  be  necessary  to  know,  for 
instance,  the  conditions  of  the  arc  prior  to  the  time  the  acceleration 
begins  and  this  information  is  not  currently  available.  Probably  so 
detailed  a  calculation  is  not  justified  considering  our  current  state 
of  experimental  knowledge.  In  future  experiments  it  would  be  desirable 
if  investigators  could  measure  the  temperature  of  the  arc  by,  for  ex¬ 
ample,  spectroscopic  means.  The  adequacy  of  various  theoretical  treat¬ 
ments  could  then  be  ascertained  and  more  sophisticated  calculations 
undertaken. 

Some  of  the  remaining  assumptions  made  in  the  model  calculation 
merit  some  brief  discussion.  First,  we  found  that  over  most  of  the  arc 
essentially  all  of  the  ions  were  doubly  ionized.  The  high  degree  of 
ionization  suggests  that  some  triple  ionization  is  probably  also  present, 
and  triple  ionization  is  not  accounted  for.  As  has  been  pointed  out, 
however,  in  a  two-  or  three-dimensional  model,  one  expects  lower  tempera¬ 
tures  and,  in  that  case,  higher  ionization  is  probably  negligible.  Ac¬ 
count  could  be  taken,  of  course,  of  triple  ionization  by  including  an 
additional  equation  in  the  hierarchy  represented  by  Eqs.  (3.53).  The 

21 

spectrum  of  trebly  ionized  copper  has  been  reported  in  the  literature  , 
and  so  the  calculation  should  be  quite  straightforward.  In  practice, 
however,  the  algebra  becomes  extremely  tedious  and  is  probably  not  justi¬ 
fied  in  view  of  the  rather  limited  improvements  anticipated  in  the  re¬ 
sults  . 


The  relative  importance  of  ordinary  heat  conduction,  neglected  in 
the  calculation,  and  radiation  can  be  estimated  by  comparing  the  co¬ 
efficient  of  thermal  conductivity  to  the  coefficient  of  radiation  heat 
conduction.  The  latter  coefficient  is  defined  in  Eq.  (3.38),  viz., 

16a  X  T3 

kr  =  — I - >  t5-4^ 


whereas  the  thermal -conductivity  coefficient  can  at  least  be  approxi¬ 
mated  by  the  expression  appropriate  for  a  Lorentz  gas.  One  has"^ 
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21.  J.M.  Schroeder  and  Th.  A.M.  Van  Kleef,  "The  Spectrum  of  Trebly 
Ionized  Copper",  Physica  49,  388  (1970) . 


55 


and  the  ratio  of  the  two  values  is 


K  R  9 

—  =  1.5  x  1(T 
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1 .23  X  1Q7T5//2 

Z  /n“ 
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Radiation  mean  free  paths,  X,  in  the  calculation  were  found  to  be  typi¬ 
cally  of  the  order  of  10  ^m,  and,  assuming  as  approximate  mean  values, 
Z=2,  T  =  5.6  X  10  deg,  n  -  10  ^  m  we  find  k^/k  ^  15.  Thus,  radi- 

ation  is  substantially  more  important  than  heat  conduction  for  the  case 

under  study  here.  The  small  value  of  X  quoted  above,  incidentally, 

serves  also  to  justify  representing  the  radiation  flux  by  Eq.  (3.38). 

As  pointed  out  in  Sec.  IIIB,  the  validity  of  that  assumption  depends  on 

the  condition  X  <<  £  . 

a 


It  was  also  suggested  in  Sec.  IIIB  that  the  condition  X  <<  £ 

a 

should  imply  that  the  properties  of  the  arc  are  nearly  independent  of 
the  temperature  at  the  projectile  surface.  This  assertion  was  used  to 
justify  employing  the  boundary  condition  in  Eq.  (3.39)  to  determine  the 
surface  temperature.  To  check  the  validity  of  the  assertion,  we  have 
performed  the  calculation  holding  the  projectile  surface  at  some  dif¬ 
ferent  temperature  and  observed  little  change  in  the  results.  In  fact, 
although  we  lowered  the  tempbrature  at  the  surface  by  more  than  an  order 
of  magnitude  from  the  value  predicted  in  the  original  calculation,  we 

observed  negligible  change  in  T  within  the  interior  of  the  arc.  There 
was  a  substantial  difference  in  the  temperature  within  a  mean  free  path 
or  so  of  the  boundary  for  the  two  calculations,  but  this  difference  is 
of  little  practical  importance. 


One  of  the  least  satisfactory  aspects  of  the  present  calculation 
is  the  reliance  on  Eqs.  (4.15)  and  (4.17)  for  the  conductivity  and  mean 
free  path  for  the  arc.  As  pointed  out,  the  derivation  of  Eq.  (4.17)  is 
based  on  the  hydrogenic  approximation  and  the  resulting  expression  is 

probably  accurate  to  within  no  better  than  an  order  of  magnitude^. 
Furthermore,  the  theory  used  to  develop  Eq.  (4,15)  is  known  to  break 
down  for  high  electron  densities  and  low  temperatures.  For  the  cases 
under  study  here,  we  are  probably  very  near,  if  not  beyond,  the  limit 

22 

of  validity  of  the  theory.  Experiments  ,  however,  have  suggested  that 
Eq.  (4.15)  may  be  more  generally  applicable  than  would  be  expected  on 

purely  theoretical  grounds^.  At  any  rate,  future  attempts  to  determine 
more  accurate  expressions  for  these  quantities,  or  to  establish  more 
firmly  the  validity  of  existing  ones,  would  be  desirable. 


22.  S.C.  Lin,  E.L.  Resler  and  A.  Kantrowitz,  "Electrical  Conductivity 
of  Highly  Ionized  Argon  Produced  by  Shock  Waves",  J.  Appl.  Phys . 
26,  95  (1955). 
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It  may  be  mentioned,  incidentally,  that  a  modification  to  Eq. 

16  23 

(4.15)  has  been  derived  *  for  the  case  in  which  conduction  takes 
place  normal  to  a  strong  magnetic  field.  The  revised  expression 
accounts  for  the  tendency  of  the  electrons  to  spiral  in  such  a  field, 
reducing  the  conductivity.  Since  the  fields  under  consideration  here 
are  quite  large,  it  might  appear  that  the  modified  conductivity  would 
be  appropriate.  It  must  be  remembered,  however,  that  the  density  is 
also  quite  large  over  most  of  the  arc  and,  thus,  the  mean  free  path 
is  small.  Therefore,  the  spiralling  effect  is  substantially  reduced. 
Quantitatively,  the  condition  which  must  be  satisfied  in  order  to  use 
the  high-field  conductivity  is  that  the  electron  gyration  radius  be 
small  compared  to  its  mean  free  path.  Our  estimates  indicate,  however, 
that  the  ratio  of  these  two  parameters  is  10,  or  greater  except  very 
near  the  trailing  edge  of  the  arc.  Therefore,  the  use  of  Eq.  (4.15), 
without  the  high-field  approximation,  is  justified.  This  point  has 

7 

been  discussed  further  by  McNab  . 

The  foregoing  discussion  makes  clear  the  need  for  some  specific 
future  studies.  First,  it  would  be  worthwhile  to  extend  the  calculations 
to  two  or  three  dimensions.  As  pointed  out  previously,  the  one-dimen¬ 
sional  model  probably  overestimates  the  temperature  of  the  arc,  and  does 
not  account  for  any  spatial  variation  of  the  flow  parameters  in  direc¬ 
tions  normal  to  the  acceleration  direction.  In  addition,  use  of  in¬ 
finitely  high  rails  overestimates  the  inductance  per  unit  length  of  the 
rails  and  is  responsible,  in  part,  for  the  relatively  high  theoretical 
value  for  the  acceleration  obtained  in  Sec.  IVC.  Second,  some  consider¬ 
ation  of  timeJdependent  effects  would  be  worthwhile,  if  only  to  better 
justify  the  steady-state  assumption.  Third,  it  would  be  of  interest  to 
carry  out  the  calculations  for  arcs  other  than  the  copper-vapor  arc 
studied  here.  From  such  investigations  we  could  determine  how  the  prop¬ 
erties  of  the  arc  material  affect  the  overall  arc  dynamics  and  which 
arcs  produce  the  most  desirable  effects.  Finally,  some  effort  should  be 
expended  to  determine  the  amount  of  heat  transferred  to  the  rails. 

When  compared  to  the  results  for  more  conventional  guns,  the  analysis 
should  provide  a  preliminary  step  in  the  study  of  gun-tube  erosion  for 
the  rail  gun. 
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APPENDIX 


The  purpose  of  this  appendix  is  to  prove  the  statement  made  in 
Sec.  Ill  that  the  acceleration  of  the  arc  and  projectile  is  independent 
of  how  the  current  is  distributed  in  the  arc  and  is  given  by 


From  Eqs.  (3.5)  and  (3.29)  we  have 


liA 
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+p* 
a  i 


J(5) 


J  J(5')d5' 
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Now  define  the  function 


u(C) 


J(5)d5 


(A.  3) 


and  differentiate  to  produce 


du 

d? 


=  -  J(C). 


(A. 4) 


Substituting  Eqs.  (A. 3)  and  (A. 4) 


into  Eq.  (A. 2),  we  find 


(A.  5) 


whence 


However,  according  to  Eq.  (3.16),  u(0) 


(A.  6) 


~  ,  whereas  u(l)  =  0,  so 
a 


(A.  7) 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR^TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1 .  BRL  Report  Number 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


